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PRECEDING  PAO!  BIXNK-NOT  FILMED 


ON  THE  RESOLVABILITY  OF  SINUSOIDS  OF 
NEARLY  EQUAL  FREQUENCY 

by 

R.V.  Webber 


ABSTRACT 

The  possibility  of  resolving  two  sinusoids  has 
been  investigated  for  the  case  where  AfAx  < 1 where 
Af  is  the  difference  in  the  frequencies  of  the  two 
sinusoids  and  Ax  is  the  time  interval  during  which 
data  may  be  taken.  The  chances  of  successful 
resolution  is  very  dependent  on  the  difference  in 
the  phase  constant  of  the  two  sinusoids.  The 
chances  are  best  when  the  sinusoids  are  in  phase 
and  are  worst  when  the  sinusoids  are  out  of  phase. 
Noise  makes  resolution  more  difficult.  The  method 
of  investigation  consisted  of  drawing  contours  of 
the  rms  residual  resulting  from  the  fitting  of  two 
sinusoids,  by  least  squares,  to  the  data.  With  no 
noise  in  the  data  the  procedure  corresponds  to  find- 
ing the  Fourier  series  which  best  fits  the  data. 

With  noise  in  the  data  the  series  so  found  is  an 
approximation  only  to  the  Fourier  series  for  the 
data. 


1.  INTRODUCTION 

This  report  deals  with  the  resolution  of  two  sinusoids  of  nearly  equal 
frequency;  that  is,  the  decision  as  to  whether  there  are  two  sinusoids  in  a 
signal  or  only  one  and  the  determination  of  the  properties  of  each  signal. 

The  bulk  of  the  report  is  devoted  to  studying  the  difficulties  with  which 
any  decision  rule  must  deal.  A possible  decision  rule  is  given  at  the  end. 
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We  first  discuss  how  the  problem  arose  and  some  of  the  restrictions 
under  which  the  problem  must  be  solved.  We  then  show  that  the  approach  is  a 
special  case  of  a general  Fourier  analysis.  The  method  of  least  squares  best 
fit  is  then  described.  Following  that  we  describe  the  surface  in  frequency 
space  on  which  a search  must  be  made  for  the  frequencies  which  give  the 
minimum  rms  residual.  Contours  of  these  surfaces  are  shown  for  the  fitting 
of  two  sinusoids  to  both  a one-sinusoid  signal  and  a two-sinusoid  signal. 

The  effect  on  this  surface  of  noise  in  the  data  is  shown  as  well  as  the  effect 
of  varying  the  signal  parameters.  The  parameters  of  the  derived  signals  are 
discussed  briefly,  and  finally,  a possible  way  is  described  for  deciding 
between  a one-sinusoid  signal  and  a two-sinusoid  signal. 

The  problem  arises  from  a consideration  of  the  resolution,  in  either 
range  or  direction,  of  two  closely  spaced  radar  targets  when  a swept  FM 
sounder  is  to  be  used.  Consider  the  range  problem  first.  With  one  target,  a 
swept  frequency  CW  sounder  generates,  after  demodulation,  a sinusoid, 

Ae i*2f t + 

, whose  frequency,  f,  is  proportional  to  the  range  of  the  target. 
The  amplitude.  A,  is  proportional  to  the  strength  of  the  echo.  The  phase 
constant,  4>,  is  unknown.  With  two  targets  the  vector  sum  of  two  sinusoids  is 
generated.  The  difference  in  the  frequencies  of  the  two  sinusoids  is  propor- 
tional to  the  difference  in  the  ranges  of  the  two  targets.  Thus,  the  problem 
of  resolving  the  ranges  of  two  closely  spaced  targets  becomes  the  problem  of 
resolving  the  frequencies  of  two  sinusoids  whose  frequencies  are  nearly 
equal.  Now  consider  the  problem  of  resolution  in  direction  of  arrival.  With 
a linear  array  of  antennas  and  a swept  frequency  CW  sounder,  a single  target 
generates,  after  demodulation,  a spatial  sinusoid  Ae^^2^fxx  + <P0)t  where  x is 
distance  along  the  array.  The  spatial  frequency,  fx,  is  given  by  fx  = (cos 
f')/'  where  0 is  the  angle  between  the  array  axis  and  the  target.  Two  targets 
generate  the  vector  sum  of  two  spatial  sinusoids  along  the  array.  The  differ- 
ence in  the  spatial  frequencies  of  the  two  sinusoids  is  proportional  to  the 
difference  in  the  cosines  of  the  angles  of  the  two  targets.  Thus  the  problem 
of  resolving  the  directions  of  two  closely  spaced  targets  becomes  the  problem 
of  resolving  two  sinusoids  whose  spatial  frequencies  are  nearly  equal.  For 
this  application,  "time"  as  used  in  this  paper  would  be  replaced  by  "distance 
along  the  array". 

Several  workers  have  considered  the  general  problem  of  resolving  signals, 
usually  in  terms  of  resolving  radar  targets.  Woodward1  and  Seibert2  discuss 
it  in  terms  of  ambiguity  functions.  Ksienski  and  McGhee3* 4 investigated 
angular  resolution  beyond  the  Rayleigh  limit.  Nilsson5  examined  the  optimum 
range  resolution  of  radar  signals.  Gething6  presented  a method  of  resolving 
the  elevation  angles  of  several  incoming  signals.  Other  relevant  papers  are 
those  by  Helstrom7,  Buck  and  Gustincic9  and  Lichtenstein  and  Young9. 

It  is  assumed  that  the  data  are  digitized  and  complex,  that  is,  each 
datum  consists  of  an  amplitude  and  a phase  (or  a real  part  and  an  imaginary 
part),  and  that  the  data  has  been  recorded  in  a limited  time  interval.  We 
shall  try  to  resolve  two  sinusoids  of  small  frequency  separation.  For 
illustration  purposes  we  shall  assume  that  the  signal  duration  is  one  second 
and  that  the  frequencies  of  the  two  sinusoids  are  separated  by  0.2  Hertz. 
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2.  SOME  COMMENTS  ON  FOURIER  ANALYSIS 


Consider  first  the  classical  discrete  Fourier  analysis.  The  discrete 
Fourier  series  which  represents  the  sampled  data,  y(t^),  Is 


N - 1 


12lSi 

N 

y(t j ) - Z zne  (1) 

n = 0 

where  tj  « jAt  Is  the  time,  N Is  the  number  of  terms  in  the  expansion  and  ^ 
Zn  - Ane^n  where  Aj,  and  4>n  are  the  amplitude  and  the  phase  constant  of  the  n 
component.  In  the  usual  discrete  Fourier  analysis,  the  complex  constants, 

Zj,,  are  selected  to  make  the  series  fit  the  data  best  In  the  sense  of  least 
squares10.  The  number  of  terms,  N,  in  the  Fourier  series  is  equal  to  the 
number  of  data  and  are  then  computed  by  the  expression  > 


Z 

n 


N - 1 

¥ y(tj)e 


(2) 


Using 

other 


this  expression  each 

Z . 
n 

The  frequency  of  the 


Z^  is  found  separately  and  independently  of  all 


th  _ . 

n component  is 


f n.  „ n 
n * L “ NAt 


(3) 


where  L in  seconds  is  the  total  time  interval  during  which  data  are  taken, 
and  the  frequency  difference,  Af,  of  any  two  successive  components  is  the  same 
as  the  resolution  and  is  given  by  the  expression 


Af  L - 1 


(4) 


Thus,  if  we  have  only  one  second  of  data  the  usual  discrete  Fourier 
analysis  will  not  resolve  components  whose  frequencies  are  separated  by  less 
than  ope  Hertz. 


Now  let  us  turn  our  attention  to  a superesolution  technique,  and  show 
its  equivalence  to  the  discrete  Fourier  analysis.  If  we  have  a Fourier 
series  in  which  N,  the  number  of  terms,  if  greater  than  the  number  of  data, 
the  formalism  can  still  be  preserved  by  adding  zeros  to  the  data  to  bring 
the  total  number  up  to  N,  and  hence  increasing  L and  the  resolution  (decreas- 
ing Af  in  Equation  4).  However,  there  are  now  too  many  unknowns  and 
restrictions  (in  this  case  the  assumption  that  the  amplitude  of  most  of  the 
terms  is  zero)  must  be  introduced  for  a solution.  Our  objective  is  to  be 
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able  to  resolve  two  sinusoids  for  which  the  product  AfL  is  as  low  as  0.2  or 
perhaps  even  0.1.  The  purpose  of  the  present  investigation  is  to  learn  the 
conditions  under  which  we  can  hope  to  reach  this  goal. 

The  approach  is  to  fit  sinusoids  to  the  data  by  the  method  of  least 
squares.  We  will  consider  only  the  cases  where  all  the  energy  is  confined  to 
a band  of  about  two  Hertz.  There  may  be  only  one  sinusoid  or  there  may  be 
two  sinusoids.  For  the  present  we  assume  that  there  is  no  noise  in  the  data. 

The  Fourier  analysis  is  now  much  simplified.  There  are  at  most  two 

non-zero  A . 

n 

We  now  find  the  two  sinusoids 


Sj  - a/<2IV  + V 


Sk  - ♦ V 


(7) 


which  gives  the  least  squares  best  fit  to  the  data.  Note  that  this  can  be 
done  without  knowing  the  value  of  N,  and  that  since  there  is  no  noise  in  the 
data,  the  rms  residue  may  be  reduced  to  a value  of  zero. 


If  there  are  two  components  they  must  be  orthogonal  in  the  interval  L. 

Using  Equation  (3)  we  see  that  their  indices  n , n must  be  such  that 

j * 


(8) 


From  Equation  (8)  we  observe  that  the  ratio  of  the  frequencies  of  any 
two  components  of  a Fourier  series  is  a rational  number. 

From  Equation  (8)  and  the  relationship  L - NAt  we  have 


“l  "k 

N “ fjdAt)  “ fk(At)  (9 

We  now  know  all  the  constants  in  Equation  (1);  that  is,  we  know  a 
Fourier  series  for  the  data.  This  series  will  be  orthogonal  in  the  interval 
NAt,  but  not,  in  general,  in  the  Interval  during  which  data  were  taken. 

With  noise  in  the  data  our  procedure  will  not  correspond  strictly  to  a 
Fourier  analysis.  The  noise  introduces  extra  components  so  that  the  Fourier 
series  must  have  more  than  two  terms.  Moreover,  the  terms  we  find  will  not 
be  Independent  of  the  other  terms  since  they  are  not  fitted  to  the  data 
throughout  the  whole  Interval  of  orthogonality. 

However,  with  a high  signal  to  noise  ratio,  the  extra  terms  will  be 
small  when  compared  with  the  signal  sinusoids  and  the  error  which  is  intro- 
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duced  by  neglecting  the  small  terms  will  be  insignificant.  This  error  will 
increase  as  the  signal  to  noise  ratio  decreases. 


3.  LEAST  SQUARES  METHOD 

Let  the  data,  y j , be  sampled  and  let  each  datum  be  a complex  number. 

For  illustration  purposes  we  assume  that  we  have  fifteen  data.  Let  Sj  and  S2 
be  the  sinusoids,  Aje^27^^  + ^01^  and  A2eK27Tf  2C  + which  we  are  to 

fit  to  the  data  by  the  method  of  least  squares.  Let 


- si  - 


(10) 


The  problem  is  to  find  Si  and  S2  such  that  Q is  a minimum.  Note  that  Q is  a 
function  of  Ai,  A2,  4>oi*  4>02»  fi,  fz  and  each  of  these  must  be  treated  as  an 
independent  variable  in  finding  the  minimum.  The  calculation  of  the  minimum 
is  done  in  two  parts. 

First  Q is  minimized  with  respect  to  the  set  (Ai , A2,  0oi»  $02)  for 
specified  values  of  fi  and  f2.  This  is  done  algebraically  by  solving  the 
linear  complex  matrix  equation 


where 


CZ  = Y 


Z = (Z  , Z ) 
1 2 


(ID 

(12) 


C = 


c c 

11  12 

c c 

21  22 


(13) 


Y = (Y  , Y ) 
1 2 


(14) 


C , - £ (X*X  ) 
j k 


m - 1 


(15) 


j m j 


m - 1 


(16) 


M - 1 


v i27Tf  At 
Xj  - e j 


(17) 
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Z1  ' V^J  <18) 

and  the  T in  Equations  (12)  and  (14)  means  "transpose".  These  equations  are 
derived  in  Appendix  A. 

Let  Q'  be  this  minimum  in  Q with  respect  to  (Ai , A?,  <£oi»  $02)  for 
specified  values  of  fj  and  f2.  It  follows  that  a minimum  in  Q'  » Q'(fi,  f2), 

or  in  with  respect  to  (fi,  fa)  coincides  with  a minimum  in  Q with  respect 

to  (Aj , A2,  <f>02>  fi»  f 2 ) • 

Thus  the  second  step  is  to  minimize  S(p  with  respect  to  (fj,  fa).  To 
find  such  a minimum  the  computer  calculates  values  of  Ap  in  the  ( f 1 , fa) 
plane  and  plots  contours.  In  this  work  contour  plots  of  v^p  were  investigated 

near  a minimum.  The  computer  program  would  first  search  for  a given  contour 

along  the  45°  diagonal.  If  found,  the  contour  would  be  traced.  If  not  found 
on  the  45°  diagonal,  a search  for  the  contour  was  made  along  the  135°  diagonal. 
The  two  diagonals  were  searched  because  within  the  circle  in  which  the  contours 
were  traced  some  contours  crossed  the  45°  diagonal  which  did  not  cross  the  135° 
diagonal  and  some  crossed  the  135°  diagonal  which  did  not  cross  the  45° 
diagonal. 

To  simulate  a signal  the  real  and  imaginary  parts  of  a one-sinusoid 
signal  or  of  the  vector  sum  of  a two-sinusoid  signal  were  generated  in  the 
computer  for  each  datum.  The  time  t = 0 was  taken  to  be  the  middle  of  the 
data  interval.  To  simulate  Gaussian  noise,  independent  and  random  Gaussian 
numbers  were  generated  such  that  the  expected  value  was  zero  and  the  expected 
rms  value  was  equal  to  the  square  root  of  the  mean  noise  power,  On.  Separate 
random  numbers  were  added  to  the  real  and  the  imaginary  parts  of  each  datum. 

The  S/N  ratio  was  found  from  the  ratio  of  the  amplitude  of  the  largest 
component,  which  was  always  one,  to  CJn,  i.e. 

S/N (dB)  = -20  log  O (19) 

1 0 n 


4.  CONTOURS  FOR  A ONE-SINUSOID  SIGNAL 

For  these  contours  the  simulated  signal  contains  one  sinusoid,  Sa,  and 
noise  when  specified.  Two  sinusoids  Si  and  S2  are  fitted  to  the  data,  and 
the  residual  errors  are  studied. 

The  region  where  contours  are  drawn  is  the  semicircle  above  the  45° 
diagonal  as  shown  in  Figure  1.  The  radius  of  the  semicircle  corresponds  to 
0.3  Hz.  This  would  correspond  in  operational  terms  to  knowing  that  wave 
energy  was  confined  to  within  0.3  Hz  of  fa,  and  looking  for  two  frequencies 
in  this  band.  No  contours  are  drawn  below  the  diagonal  since  the  diagram  is 
symmetric  in  fi  and  f2.  The  frequency  of  the  signal,  fa,  is  represented  by 
the  two  orthogonal  dashed  lines.  This  would  be  the  oblective  of  our  search 
in  an  operational  situation. 
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Figure  I.  Diagram  showing  semi-circle  in  which  contours  were  traced  for  a one-sinusoid  signal 


The  numbers  on  the  contours  in  the  remaining  figures  are  the  height  of 
the  contour  in  dB  with  respect  to  the  amplitude  of  the  signal.  The  two 
orthogonal  straight  lines  are  fi  » fa,  £2  * fa.  The  contours  for  which  a 
search  was  made  were,  in  dB,  -54,  -46,  -40,  -30,  -26,  -20,  -14,  -10  and  -8. 
For  all  plots  for  a one-sinusoid  signal  the  amplitude  of  the  signal  was 
unity  and  the  phase  constant  was  zero.  Other  values  of  the  phase  constant 
yielded  identical  plots. 


4.1  NOISE  FREE  DATA 

In  Figure  2 contours  are  plotted  for  a pure  one-sinusoid  signal  which 
contains  no  noise.  The  height,  v^Q*",  of  the  contours  with  respect  to  the  signal 
amplitude  are  given  in  dB  by  the  expression 

dB  - 20  log  ^Qr 
1 0 

If  either  fi  or  £2  equals  £a  the  rms-residual  is  zero.  If  the  rms-residual 

is  near  zero  then  one,  or  both,  of  fi  or  £2  will  be  near  fa.  However,  we 

' 


j 
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cannot  without  additional  checking  use  this  characteristic  to  distinguish 
between  a one-sinusoid  signal  and  a two-sinusoid  signal  as  can  be  seen  by 
considering  representative  points.  At  points  A,  B and  C of  Figure  2 the 
amplitudes  of  the  derived  signals  are  as  sketched  in  Figure  3,  A,  B and  C 
respectively.  These  points  are  near  the  zero  rms-residual  lines  and  are 
representative  of  the  end  points  of  a search.  At  points  A and  C of  Figure  2 
the  amplitude  which  corresponds  to  the  fitted  frequency  which  is  closest  to 
fa  is  large  while  the  other  is  small.  From  this  one  might  conclude  that  the 
signal  was  made  up  of  a main  component  and  possibly  another  small  component. 
At  point  B however,  the  two  fitted  sinusoids  have  equal  amplitudes  from  which 
one  might  conclude  that  the  signal  is  made  up  of  two  sinusoids  of  equal 
amplitudes  - which  would  be  incorrect. 


4.2  EFFECT  ON  NOISE  IN  THE  DATA 

The  effect  on  the  rms-residual  contours  of  adding  noise  to  the  data  is 
shown  in  Figures  4 to  6 where  the  S/N  ratios  are  50  dB,  40  dB  and  30  dB 
respectively.  As  the  S/N  ratio  decreases  the  lower  contours  are  affected  by 
noise  but  for  the  most  part  they  hold  their  shape  and  position  until  they 
disappear  entirely. 

Of  the  contours  which  were  made  to  disappear  by  the  noise,  the  highest 
was  equal  to,  or  about  equal  to,  the  rms  amplitude  of  the  noise  itself.  The 
rms-residual  could  be  reduced  only  to  within  about  6 dB  of  the  rms  amplitude 
of  the  noise. 

Figures  7 and  8 give  some  idea  of  how  much  a contour  varies  in  shape 
and  position  with  different  samples  of  noise  in  the  data.  For  these  figures 
the  same  contour  was  drawn  10  times.  Each  time  the  signal  was  the  same  but 
with  a different  set  of  noise  samples.  For  Figure  7 the  S/N  ratio  was  50 
dB.  For  Figure  8 the  S/N  ratio  was  40  dB.  All  contours  were  confined  to 
the  shaded  areas  of  the  figures. 


4.3  EFFECT  OF  VARYING  SIGNAL  FREQUENCY 

Figure  9,  10  and  11  are  for  a signal  with  no  noise  but  with  the  signal 
frequency  varying  linearly  in  time  by  0.01  Hz  sec  1 , 0.04  Hz  sec  1 and  0.08 
Hz  sec  1 respectively.  The  general  effect  is  the  same  as  having  noise  in  the 
data;  the  higher  variation  in  frequency  corresponding  to  the  lower  signal-to- 
noise  ratio. 


5.  CONTOURS  FOR  A TWO-SINUSOID  SIGNAL 

For  the  two-sinusoid  signal  the  area  for  which  contours  were  drawn  is  a 
segment  of  a circle  above  the  45°  diagonal  centred  on  the  two  frequencies, 
fa,  fg  in  the  simulated  signal.  This  is  illustrated  in  Figure  12.  Presum- 
ably, in  an  operational  situation,  a semicircle,  centred  on  the  centre  of  the 
frequency  band  would  be  used.  The  same  contours  were  drawn  as  for  the 

single-sinusoid  signal.  Unless  otherwise  stated  A - A,,  = l . 

u p 


# 


Figure  2.  Rms-residual  contours  for  fitting  two  sinusoids  to  a one-sinusoid  signal.  No  noise  in  signal. 


f. 


fi 

(A) 


fa 


fl 

(C) 


Figure  3.  Diagram  showing  possible  relationships  of  f , and  f2  to  f The  rms-residual  is  the  same  in  all 
three  cases.  The  lengths  ot  the  arrows  represent  amplitude.  Cases  A,  B and  C correspond  to  the  points 

A,  B and  C of  Figure  2. 
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Figure  12.  Diagram  showing  semi<irde  in  which  contours  were  traced  for  a two-sinusoid  signal. 

The  orthogonal  axes  in  the  following  diagrams  are  fj  * fa,  f2  = fg. 

The  fitting  of  a single  sinusoid  to  the  data  will  give  a point  on  the  line 
fi  * f2,  of  Figure  12.  To  distinguish  between  a one-sinusoid  signal  and  a 
two-sinusoid  signal  we  must  be  able  to  distinguish  between  a point  on  the 
diagonal  and  a point  off  the  diagonal.  In  an  operational  situation,  contours 
would  be  plotted  with  the  objective  of  finding  the  origin  in  Figure  13  and 
onward.  More  contours  indicate  steeper  slopes  of  the  surface  VQ*  in  (fi,  f2) 
space  and  hence  a greater  probability  of  finding  the  minimum. 

5.1  NOISE  FREE  DATA  - EFFECT  ON  CHANGING  FREQUENCY  SEPARATION  AND  PHASE 

Figure  13  shows  the  contours  when  fa  - 8.968  and  fg  - fa  “ 0.4  Hz  for 
4>g  - <>a  “ 90°  and  180°.  The  shape  of  the  contours  is  very  dependent  on 

the  phase  difference,  <f>g  - <pa.  It  would  appear  that  the  chances  of  resolving 
two  sinusoids  would  be  best  when  the  two  signal  sinusoids  are  in  phase  since 
it  is  for  this  case  that  the  point  representing  the  signal  is  separated  from 
the  line  fi  - f2  (the  flat  side  of  the  circle)  by  the  largest  number  of 
contours.  At  90°  out  of  phase  the  chances  would  be  less  since  some  of  the 
contours  are  open  towards  the  45s  diagonal.  At  180°  out  of  phase  even  more 
of  the  contours  are  open.  The  chances  of  resolving  two  sinusoids  is  poorest 
if  they  are  180°  out  of  phase.  Figure  14  is  a similar  plot  for  the  same 
value  of  f2  but  for  fg  - fa  ■ 0.2  Hz.  The  contours  are  larger  and  the  sur- 
face is  less  steep.  The  point  which  represents  the  signal  is  less  well 
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separated  from  the  diagonal;  this  makes  resolution  more  difficult.  For 
<{>g  - 4>a  = 180°  even  the  -54  dB  contour  is  open.  One  would  have  to  go  to 
lower  contours  to  find  one  that  is  closed  - perhaps  to  -60  dB  or  -70  dB. 

Figure  15  is  again  a similar  plot  for  the  same  value  of  fa  but  for 
fg  - fa  = 0.1  Hz.  Again  the  contours  are  larger  and  the  surface  flatter. 

For  = 180°  one  would  have  to  go  to  still  lower  contours  to  find  one 

that  is  closed. 

There  would  be  two  possible  ways  to  use  these  contours  to  aid  in  finding 
the  approximate  position  of  the  center  of  the  axis. 

1.  To  search  the  surface  until  the  rms  residual  is  below  a certain 
contour  level.  One  would  then  know  that  the  point  found  would  be 
inside  that  contour.  Clearly  this  would  only  prove  the  existence 
of  two  sinusoids  if  the  point  found  were  on  a closed  contour,  and 
would  probably  not  give  their  properties  very  accurately. 

2.  Draw  rms-residual  contours.  One  could  then  choose  the  "center"  of 
the  contour.  This  would  give  greater  accuracy  than  method  (1). 

Also  contours  which  were  not  closed  could  be  used.  It  would,  how- 
ever, be  a more  lengthy  procedure  than  that  of  method  (1). 

Of  these  two  methods,  either  one  could  be  used  if  the  two  signal 
sinusoids  are  in  phase  or  even  if  they  are  90°  out  of  phase.  With  the  two 
signal  sinusoids  180°  out  of  phase,  however,  method  (1)  would  place  such 
demands  on  the  S/N  ratio  that  it  would  be  useless  in  most  cases  for  Af  0.2. 
There  is  therefore,  little  hope  of  resolving  two  sinusoids  with  AfAt  0.2  if 
they  are  almost  180°  out  of  phase.  And  this  is  with  no  noise  in  the  data. 

It  was  found  that  the  contours  depended  only  on  the  difference  between 
the  signal  frequencies  and  not  on  the  absolute  frequencies. 


Injure  13.  Rms-residual  contours  for  fitting  two  sinusoids  to  a two-sinusoid  signal.  Ip  - fa  = 0.4  H/;  tor 
A,  B and  C yg  - - (f,  9(f  and  180°  respectively.  No  noise.  The  lowest  < uniuws  are  at  -54  dB. 
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figure  14.  Rms-residual  contours  tor  fitting  two  snusotds  to  a two-sinusoid  signal,  Ip-  lQ=  0.2  Hi;  for 
A,  H and  < *Pp  - <Pa  = 0°,  90°  and  180°  respectively.  No  noise.  The  lowest  contours  are  at  -5 4 dB. 


Figure  15.  Rms-residual  contours  for  fitting  two  sinusoids  to  a two-sinusoid  signal,  fp-  tQ  = 0. 1 Hi;  for 
A,  B and  C g>p  - <Pa  = <f,  90°  and  1 80°  respectively.  No  noise.  The  lowest  contours  in  A and  C are 
-54  dB;  the  lowest  in  B is  -46  dB.  Because  of  computer  program  failure  the  -54  dB  contour  was 

missed  in  B. 


5.2  EFFECT  OF  NOISE  IN  THE  DATA 

Figures  16,  17,  18  and  19  show  what  happens  to  the  contours  for  which 
fg  - fa  * 0.2  Hz  when  the  S/N  ratio  is  50  dB,  40  dB,  30  dB  and  20  dB 
respectively.  Again  as  the  noise  Increases  the  lower  contours  disappear  and 
the  rms  residual  can  be  decreased  only  to  within  6 dB  of  the  rms  amplitude  of 
the  noise.  At  the  same  time  the  higher  contours  contract  a little  and  some 
shifting  of  the  contours  is  evident  at  the  lower  S/N  ratios  but,  for  the  most 
part,  these  contours  are  not  changed  very  much. 

In  the  figures  thus  far  for  two  sinusoid  signals  each  diagram  corres- 
ponds to  one  set  of  noise  samples,  one  sample  for  each  datum.  To  find  the 
variation  of  the  contours  as  the  noise  samples  change  the  same  contour  was 
traced  ten  times;  each  time  the  S/N  ratio  was  the  same  but  the  noise  samples 
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were  different.  Figures  20  and  21  show  the  results  for  S/N  ratios  of  AO  dB 
and  30  dB,  respectively.  At  a high  S/N  ratio  the  location  of  the  contours 
does  not  move  much;  as  the  S/N  ratio  decreases  the  curves  move  around  more. 

We  may  note,  however,  that  even  at  S/N  ” 30  dB  when  the  two  signal  sinusoids 
are  in  phase  or  90°  out  of  phase,  the  "centers"  of  all  the  contours  are  within 
0.1  Hz  of  the  point  which  corresponds  to  the  signal. 


5.3  EFFECT  OF  VARYING  FREQUENCY 

Figure  22  shows  contours  for  no  noise  in  the  data  but  for  both  of  the 
signal  frequencies  varying  by  0.01  Hz  sec  *.  The  lower  contours  have 
disappeared.  The  effect  is  very  similar  to  adding  noise  to  the  data. 


Figure  16.  Rms-residual  contours  for  fitting  two  sinusoids  to  a two-sinusoid  signal,  fa  - fQ  = 0.2  Hi;  for 
A,  B and  C <pp  - <pQ  = 0°,  90°  and  180°  respectively.  5//V  = 50  dB.  Lowest  contours  are  at  -40  dB. 


Figure  1 7.  Rms-residual  contours  for  fitting  two  sinusoids  to  a two-sinusoid  signal,  fp-  fa*  0.2  Hz;  for 
A,  B and  C <i>p  - <pQ  = 0°,  90°  and  18(f  respectively.  S/N  = 40  dB.  Lowest  contours  are  at  -30  dB. 
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C 

Figure  18.  Rms-residuai  contours  for  fitting  two  sinusoids  to  a two-sinusoid  signal,  fa  - tQ  - 0.2  Hz;  for 
A,  R and  C <Pq-  <Pa  = CF , 90  and  180°  respectively.  S/N  = 30  dB.  Lowest  contours  are  at  - 26  dB 


A 


Figure  19.  Rms-residuai  contours  for  fitting  two  sinusoids  to  a two-sinusoid  signal,  fa  - f - 0.2  Hz;  for 
A,  B and  C 4>q-  Qa  = 0°,  90°  and  180°  respectively.  S/N  = 20  dB.  Lowest  contours  A and  B are  at  -14  dB. 

No  contours  found  in  C. 


figure  20.  Rms-residuai  contours  for  fitting  two  sinusoids  to  a two-sinusoid  signal,  fp-  fQ  = 0.2  Hz;  for 
A,  B and  C <t>p  - 4a  = (f,  90°  and  18(f  respectively.  Cross-hatched  area  shows  the  extent  to  which  a 
single  contour,  -26  dB,  moved  around.  S/N  = 40  dB. 
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figure  21.  Rms-residual  c on  tours  for  fitting  two  sinusoids  to  a two-sinusoid  signal,  Ip-  ta  = 0.2  Hz;  for 
A,  H and  C Qp-  0a~  0°,  90°  and  180"  respectively.  S/N  = 30  dB.  The  cross-hatched  area  shows  the 
extent  to  which  a single  contour,  -26  dB,  moved  around. 
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Figure  22.  Rms-residual  contours  for  fitting  two  sinusoids  to  a two-sinusoid  signal,  fn  - fa  = 0.2  Hz;  for 
A,  B and  C <Pp  - <Pa  = 0°,  90°  and  I8(f  respectively.  No  noise.  Both  fQ  and  fp  vary  by  0.01  Hz  sec"1. 
Lowest  contour  in  A is  -30  dB;  in  B -40  dB  and  in  C -46  dB. 


5.4  EFFECT  OF  VARYING  AMPLITUDE  RATIO 

Figures  23  and  24  show  contours  for  no  noise  in  the  data  and  for 
- 1.0  but  for  Ag  ■ 0.5  and  0.25  respectively.  The  contours  are  stretched 
along  the  axis  which  correspond  to  the  signal  with  the  reduced  amplitude. 

The  lower  the  amplitude  the  greater  is  the  stretching.  Stretching  corresponds 
to  reducing  the  slope  in  the  direction  of  the  stretch.  This  makes  resolution 
more  difficult  in  that  direction. 

1 

5.5  ANALYSIS  OF  TWO  SUCCESSIVE  REALIZATIONS  AT  THE  SAME  TIME 

Two  time-separated  sets  of  measurement  of  a two  sinusoid  signal  may  be 
analyzed  together  in  such  a way  that  the  rms  residual  of  the  two  realizations 
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hiyure  23.  Rms-residuul  contours  for  fitting  two  sinusoids  to  a two-sinusoid  signal,  fg-  fQ  = 0.2  Hr;  for 
A,  ti  and  C 0g-  *>a  = 0°,  90°  and  18Cf  respectively.  No  noise.  A&  = 1.0;  Ag=  0.5.  * Lowest  contours 

are  at  -54  dB. 


figure  24.  Rms-residual  contours  for  fitting  two  sinusoids  to  a twc-sinusoid  signal,  fg  - f = 0.2  Hz;  for 
A,  B and  C <Pg-  9a=  0°,  90°  and  I8(f  respectively.  No  noise.  Aa  = 1.0;  An-  0.25*  Lowest  contours 

are  at  -54  dB. 


is  a minimum.  The  difference  between  the  phase  constants  of  the  two  sinusoids 
in  the  first  realization  will,  in  general,  be  different  from  that  in  the 
second.  The  result  is  identical  to  the  result  for  a single  realization  whose 
difference  in  the  phase  constant  is  intermediate  between  the  two  which  are 
processed.  Since  the  maximum  difference  in  the  phase  constants  is  180°  the 
difference  in  the  phase  constants  of  the  equivalent  single  realization  will, 
in  general,  be  less  than  180°. 

For  example,  two  realizations  were  processed  together,  one  for  which 
♦g  - 4>a  - 0°  the  other  for  which  “ ♦a  ■ 90°.  The  contours  which  were 
obtained  were  Identical  to  those  which  were  obtained  by  processing  a single 
realization  for  which  4>g  - 4>a  " 60®. 

Equations  which  relate  the  equivalent  realization  to  the  two  realiza- 
tions which  are  processed  are  derived  in  Appendix  B. 
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By  processing  two  realizations  together  we  could  hope  to  improve  the 
resolvability  of  two  sinusoids  for  it  is  much  less  likely  that  two  realiza- 
tions will  both  be  near  to  180°  out  of  phase  than  it  is  that  one  will  be.  The 
equivalent  single  realization  is,  therefore,  less  likely  to  be  near  to  180° 
out  of  phase  than  is  a single  realization. 

A disadvantage  of  processing  two  realizations  together  is  that  two  data 
intervals  are  required.  Thus  the  time  is  doubled  and,  for  a given  Af,  AfAx 
is  doubled  where  Ax  is  the  total  time  during  which  data  are  taken.  This  tends 
to  nullify  any  advantage  that  might  be  gained  since  one  of  our  main  objectives 
is  to  keep  AfAx  small. 


5.6  THE  AMPLITUDES  OF  THE  FITTED  SINUSOIDS 


The  derived  amplitudes  of  the  fitted  sinusoids  vary  in  a complex  manner. 
They  depend  on  (1)  the  relation  of  the  frequencies  of  the  fitted  sinusoids  to 
the  frequencies  of  the  signal  sinusoids,  (2)  the  difference  between  the  phase 
constants  of  the  signal  sinusoids,  (3)  the  height  of  the  contour  being 
followed  and  (A)  the  amount  of  noise  in  the  data. 


Consider  the  case  for  which  the  amplitudes  of  the  sinusoids  in  the 
simulated  signal  are  equal  and  with  no  noise.  For  this  case  it  was  found 
that  the  amplitudes  of  the  fitted  sinusoids  are  practically  always  unequal. 
Sometimes  the  one  and  sometimes  the  other  may  be  larger.  Both  may  be  smaller 
than  the  amplitudes  of  the  signal  sinusoids  or  both  may  be  larger,  by  several 
hundred  per  cent,  or  one  may  be  smaller  and  the  other  larger.  However,  they 
are  always  about  equal  on  the  135  - 315°  diagonal.  On  the  135°  radius  this 
corresponds  to  the  relation  between  (f  , fQ,  f , f ) as  sketched  in  Figure  25a. 
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Figure  25. 


Where  fi  - fa  is  negative  and  f2  - fg  is  positive  but  both  are  of  equal 
magnitude.  On  the  315°  radius  this  corresponds  to  the  relation  shown  in 
Figure  25b  where  f\  - fa  Is  positive  and  f2  - fg  is  negative  and,  again, 
both  are  of  equal  magnitude.  Of  course,  fewer  contours  exist  on  the  315° 
radius  than  on  the  135a  radius. 


Noise  in  the  data  caused  variations  in  the  amplitudes  of  the  fitted 
sinusoids  but  the  general  pattern  of  the  results  was  unchanged. 
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5.7  THE  PHASES  OF  THE  FITTED  SINUSOIDS 

5.7.1  Noiseless  Data 

When  the  two  signal  sinusoids  are  in  phase  the  two  fitted  sinusoids  are 
always  in  phase  with  the  signal  sinusoids  for  the  lower  contours  which  are  the 
important  ones  since  a search  presumably  ends  up  on  these.  For  higher  contours 
the  two  fitted  sinusoids  may  be  in  phase  with  the  signal  sinusoids  or  one  may 
be  in  phase  and  the  other  out  of  phase.  The  shift  from  in  phase  to  out  of 
phase  always  occurs  off  the  contour  diagram.  It  is  suggested  that  these  higher 
contours  have  more  than  one  branch  and  that  for  a given  branch  the  fitted 
sinusoids  are  always  in  phase  or  are  always  out  of  phase. 

When  the  signal  sinusoids  are  90°  out  of  phase,  the  phase  difference 
between  the  fitted  sinusoid  varies  from  47°  to  170°  for  the  higher  contours 
where,  again,  there  appears  to  be  more  than  one  branch  to  a contour.  For  the 
lower  contours  their  phase  difference  is  closer  to  90°. 

When  the  signal  sinusoids  are  180°  out  of  phase  the  fitted  sinusoids 
are  always  in  phase  with  them. 

5.7.2  Noisy  Data 

Putting  noise  in  the  data,  such  that  the  S/N  ratio  was  as  low  as  30  dB, 
had  little  effect  on  the  phases  of  the  fitted  sinusoids  when  the  signal 
sinusoids  were  either  in  phase  or  180°  out  of  phase.  When  the  signal  sinu- 
soids were  90°  out  of  phase,  noise  had  some  effect  on  the  phases  of  the  fitted 
sinusoids  but  the  general  pattern  of  the  results  was  unchanged. 


6.  DECISION  RULE 

As  discussed  in  Section  4.1  we  cannot  use  the  minimum  rms  residual 
alone  to  decide  between  a one-sinusoid  signal  and  a two-sinusoid  signal.  A 
possible  way  of  making  this  decision,  on  the  basis  of  the  work  reported  here, 
would  be  to  draw  the  rms  residual  contour  map  and  to  compare  it  by  eye  with  a 
typical  one-sinusoid  map  and  a typical  two-sinusoid  map  as  given  in  this 
paper.  If  the  pattern  appears  to  be  similar  to  a two-sinusoid  map  then  the 
"center"  of  the  contours  can  be  estimated,  again  by  eye.  The  vertical,  or 
horizontal,  distance  of  this  "center"  from  the  45°  diagonal  will  then  give  an 
estimate  of  the  frequency  separation  between  the  two  components  of  the  signal. 


7.  SUMMARY 

The  possibility  of  resolving  two  sinusoids  for  which  AfAt  <<  1 has  been 
investigated.  Contours  of  constant  rms-residual  were  generated  by  fitting, 
using  the  method  of  least  squares,  of  two  sinusoids  to  a one-sinusoid  signal 
as  well  as  to  a two-sinusoid  signal.  From  these  contours  It  was  concluded 
that  the  chances  of  the  successful  resolution  of  two  sinusoids  depends 
markedly  on  the  difference  of  their  phase  constants.  It  should  be  easiest 
if  they  have  the  same  phase  constants  and  most  difficult  if  their  phase  differ 
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by  180°.  Noise  eliminates  the  lower  contours  and  hence  decreases  the 
precision.  The  rms-residual  can  be  reduced  nearly  to  the  rms  noise  amplitude 
but  no  further.  The  effect  of  changing  the  ratio  of  the  amplitudes  of  the 
signal  sinusoids  and  of  changing  their  frequencies  was  investigated.  With  no 
noise  in  the  data  our  procedure  corresponds  to  the  finding  of  a Fourier 
series  which  fits  the  data.  With  noise  in  the  data,  however,  the  results 
correspond  only  approximately  to  the  Fourier  series  for  the  data. 

A possible  decision  rule  would  be  to  draw  contour  maps  of  the  rms- 
residuals  and  to  compare  them  with  typical  maps  of  a one-sinusoid  signal  and 
a two-siu', soid  signal. 
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APPENDIX  A 


Derivation  of  Equations  (11)  to  (18)  of  the  Text 


„ , . A i ( 2tt  f t +4 

S (t  ) - A e j m o j 

J m j 


t - (m  - l)At;  (m  **  1,  N) 

m 


V 12TTf  At 
Xj  - e j 


id» 

Z = A e 
j T 


where  Aj  , O0j  and  fj  are  the  amplitude,  the  phase  constant  and  the  frequency, 

respectively,  of  sinusoid  Sj(tm),  tm  is  the  time  of  datum  m(m  = 1,  N) , At 

Is  the  time  interval  between  successive,  equally-spaced  data  and  N is  the 
number  of  data  which  were  taken. 


From  equations  (Al)  to  (A4)  we  have 


W ' Vj 


m - 1 


Now  let 


P - £ |Y  - S (t  ) - S (t  ) 

m . 1 m 1 m 2 m 


where  y is  the  datum  which  was  taken  at  time  t . 

m 

We  must  find  the  value  of  Z and  Z which  make  P a minimum. 

l 2 


£ (Y  - S (t  ) - S (t  )][Y*  - S* (t  ) - S*(t  )] 
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P - £ [V  -ZXm-1-ZX  m_1]  (Y*  - Z*(X*)m  1 - Z*(X*)m  " X] 
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Let  Z - Z + C 5 (A9) 

1 1,0  11 

where  Zj  0 is  the  value  of  Zj  which  makes  P a minimum,  C i is  a small,  real 
number  and  £i  is  an  arbitrary  complex  number  which  is  not  a function  of  Ci  • 
Thus  at  the  minimum  value  of  P we  have 


- 0 

= 0 


Now 


P = 


= 1 


- z xr 

1 t 0 


- 1 


m - 1 


- Z2X2 


m - 1 ] 


(A10) 


x [Y*  - Z*  (X*)®  ‘ 1 -C!C*(X*)m  " 1 - Z2(X*)m  “ 1]}  (All) 

m i , o l li  2 

So 

3P 
*1 


L. 


UlX?  ' 1 [ Y* 


m 


Z* 
l » o 


(X*) 

1 


m - 1 


CjC*(x*) 

i i 


m - 1 


z*(x*)m  ' 1] 

2 2 


+ £*(X*)m  “ X[Y  - ZioXi®  “ 1 -CiCiXim  1 - Z2X2m  ~ X]}  (A12) 

l l ® 


2RU*  (X*)m  “ X[Y  - Z 

1 m = 1 1 m »• 


Xim  ■ 1 - CjCiXi®  ■ 1 - z2x2m  “ x]} 


(A13) 


and 


(£-)  = 2r{ f; i £ <x*)m  “ Vm  ~ z 

V *'Cl  - 0 m-1  1 m 1,0 

From  equations  (A10)  and  (A14)  we  have 


m “ 1 - z2x2m  ' 1]} 


(A14) 


II 

2R(C*  2 <X*>m  “ X[Y  - Z Xim  ‘ 1 - Z2X2m  ~ X]}  - 0 

i ,1  m l » o 

m 53  i 


(A15) 
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This  latter  equation  is  true  for  any  value  of  £,.  The  only  way  that 
this  can  be  so  is  for  the  following  equation  to  be  true. 


N 

E{(x*)m  1 y - z x ” ' 1 - nm  ' ' o 

1 m 


or 


N 

(x*x2)m  1 - £ Y (x*)m  ” 1 

1 m - 1 m 1 


(A16) 


(A17) 


Equation  (A17)  is  true  for  any  value  of  Zz. 

Now  let  Z2  » Z +C2C2 

2,0 


(A18) 


and  carry  through  the  same  procedure  again.  Equations  which  are  similar  to 
equations  (A9)  to  (A17)  will  be  obtained.  The  one  corresponding  to  equation 
(A17)  will  be 


(A19) 


z,  £ (x*Xl)m  “ 1 + z £ (x*x2)m  " 1 = £y(X*)"‘  " l. 

m = 1 m = 1 

Equation  (A19)  is  true  for  any  value  of  Zj . 

Now  choose  Z2  » Z2>0  in  equation  (A17)  and  choose  Z\  = Z^  in  equation 
(A18),  and  we  have  the  two  simultaneous  equations 


Z C11+Z  Ci2  = Yi 
l.o  2,0 


z c21  + z c2  2 = y2 

1.0  2.0 


(A20) 


where 


v • ' 1 


(A21) 


Y - £ Y (X*)“  " 1 

* m - 1 m * 


(A22) 


Solving  equations  (A20)  for  Z\  0 and  Z2>0  will  then  give  us  the  values 
of  Zi  and  Z2  which  makes  P a minimum. 
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Equations  (A20) , (A21)  and  (A22)  are  the  same  as  equations  (11)  to  (18) 
of  the  text  except  that  in  the  text  they  have  been  written  in  matrix  form  and, 
for  convenience,  the  zeros  have  been  dropped  from  the  subscripts  of  Z 

j . o 
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APPENDIX  B 


Derivation  of  Equations  for  Analyzing  Two  Realizations  Together 


Let  an  odd  number  of  data  be  taken  at  equal  Intervals  of  time,  At.  Let 
the  middle  datum  be  at  t - 0 so  that 


t * mAt 

m - -M,  - (M  + 1),  ...  -1,  0,  1,  ...  M 

(Bl) 

Let  the 

signal  be  made  up  of  the  two  sinusoids 

i(2irf  mAt  + <p  ) 
S = A e a a 

a,m  a 

(B2) 

S„  - A..1<2,,£B,*iC  * V 

p,m  p 

(B3) 

Let  an 

assumed  test  signal  be  made  up  of  the  two  sinusoids 

S - A,e1(2lTfimAt  + 

1 1 

(BA) 

s = A i(2irf2mAt  + ^ 

2 , m 

(B5) 

Let 

S - (S^  + SQ  - S - S ) 

m a, m p,m  i,m  2,m 

(B6) 

then 

| S 1 2 - 1 S |2  + |s.  |2  + 2R{S  S*  } 

1 m a.m1  B,m'  a,m  B,m 

+|S  I2  + |S  I2  + 2R|S  S*  | 

1 2 »®  1 .ro  2 ,m 


- 2R{S  S*  + S S*  + SQ  S*  + SQ  S*  1 
a,m  i ,m  cx,m  2,m  p,m  i ,m  p,m  2 , tn ; 

where  the  superscript  * signifies  the  complex  conjugate  and  R(x)  means  "real 
part  of  x". 


Let 


°Jk»  * *m<tJ ' 


(B7) 

(B8) 
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then 


S_ | 2 « + Ag  + Aj  + A2  + 2R{B„BgG„fl_  + BiB*G 


v UaU  n 1 U1  U U 

a 6 a8m  2 12m 


(B9) 


-BB*G  - B B*G  - BqB*Gq  - BoB*G0  } 
a 1 ajm  a 2 a2m  p 1 Pim  p 2 P2m 


Now 


T!  ei9m  -1+2  cos  m0 

m *■  -M 


M 

z: 

m ■ 1 


(BIO) 


Let 


H 

(>  " b?-M  Gjkm 


M 


1 + 2 co82ir(f  - f 

m - 1 J 


)At 


(Bll) 


We  note  that  G is  a real  number.  We  may  now  write 
J * 


|Sm|2  - (2M  + 1)  (A2  + A2  + A?  + Ai)  + 2 G^R^B*)  + G12R(BiB*) 


m = -M 


- G ^(B  B*)  - G R(B  B*)  - Gft  R(BftB*)  - G R(B  B*) 
oti  a 1 a2  a 2 Pi  p 1 p2  p 2 


or 


M 

E Isj’  - <2M  + 1)  (*>  + *|  + *i  + *!>  + 2(GaBAa*Bco.(4,[i  - »B) 

1 = -M 

+ Gj  2Aj A2cos (01  - 4>2) 


‘ GaiAd\c°8<*a  ” ♦»>  “ Ga2AaA2C08<<t,a  “ *2>  " G6lV‘C08(4’8  ” *l) 
- Gp  A(  A2cos(<t>g  - 4.2)}  (B12) 


Now  we  wish  to  find  the  realization,  r,  which  Is  equivalent  to  the  mean  of 
the  two  realizations  p and  q,  that  is 
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M 


M 


£j*j,1  ■ 2J,Jl  + 


(B13) 


i m - -M  “ / r \m  “ -M  “ / p \m  * -M 

From  equations  (B12)  and  (B13)  we  get 


{ (A  2 + A 2 - 2A  2)  + (Afl2  + A2  - 2A  2) 
' ap  aq  ar  Bp  Bq  Br 


+ (A  2 + A 2 - 2A  2)  + (A  2 + A 2 - 2A  2)}(2M  + 1)} 

ip  iq  ir  2P  2q  2r 


+ 2Gc16{AapABpC08(%  ' V + %ABqC0 

+ 2G  {A  A cos ((f)  - $ ) + A A cos(0 

12  IP  2 P ip  2P  iq  2q 

- 2G  {A  A cos($  - 0 ) + A A cc 

oti  ap  ip  ap  ,p  aq  jq 

- 2G  {A  A cos(0  - 0 ) + A A cc 

a2  ap  2 P ap  2P  aq  2q 

- 2GC  (A0  A cos(0fi  - 0 ) + Aq  A cc 

Bi  Bp  ip  6p  ip  Pq  iq 

- 2Gr  {Ar  A cos (0R  - 0 ) + Aft  A cc 

P2  BP  2 P t5p  2P  pq  2q  PH  2H  P*-  2 I P*-  2 

- 0 (B14) 

The  G * 8 are  functions  of  frequency  only.  Since  equation  (B14)  must  hold 
for  all  values  of  f , fQ  , f and  f each  expression  in  braces,  { },  must  be 
zero  and  we  have  011  1 1 2 

2(A  2 + A 2 + A 2 + A 2)  - A 2 + A 2 + A 2 + A 2 + A 2 + A 2 + A 2 + A 2 
ar  Br  tr  2r  ap  Bp  aq  Bq  ip  2p  iq  2q 

(B15) 

2A  A„  cos (0  - 0O  ) ■ A Aq  cos(0  - 0o  ) + A Aq  cos(0  - 0Q  ) (B16) 

ar  Br  v ar  TBr  ap  Bp  ap  TBp  aq  Bq  aq  Bq 

2A  A cos(0  - 0 ) ■ A A cos (0  - 0 ) + A A cos(0  - 0 ) (B17) 

ir  2r  ir  2r  ip  2p  ip  2P  iq  2‘1  iq  2q 

2A  A cos ( 0 - 0 ) ■ A A cos (0  “0  ) + A A cos(0  “0  ) 

ar  ir  v ar  »r'  ap  !p  -wap  ip'  aq  ,q  VTaq  yiq' 


'aq  " 

v 

< 

U 

<a 

CM 

1 

o COS (0 
Br  ar 

• v> 

« ) 

- 2A  A 

COS  (0 

• ♦ )} 

’x  q 

2q 

ir 

2r  ir 

2r 

0 ) 

- 2A  A 

COS  (0 

- 0 )} 

aq  " 

iq 

ar 

i r ar 

ir 

0 ) 

- 2A  A 

cos(0 

• 0 )) 

'aq  “ 

2q 

ar 

2r  ar 

2r 

% " 

0 „> 
iq 

- 2A3fA 

ircos(V  • 

• v» 

’6q  " 

0 n) 
2q 

- 2AgrA 

2rcos(V  ■ 

■♦«» 

(B18) 
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2AarA,rcos(V  “ V " AapA2pCOS  <%>  “ V + AaqA2qc°8<*aq  " V (Bl9) 

2ASrA1rC08(V  - V " V, pCOS1 % “ V + ABq\qcos(4>Bq  * *,q>  (B20) 

2Afl  A cos($fi  - $ ) - A0  A cos(0c  - <J>  ) + AQ  A cos(<)>  - 0 ) (B21) 

Br  2 r Br  2r  Bp  2p  Bp  2p  Bq  2q  Bq  2q 

Now  In  the  example  mentioned  in  section  V5  we  assign* d the  values. 


A - A„  » A * A„  - A » A„  - 1 
ap  Bp  aq  Bq  ar  Br 


(B22) 


In  which  case,  by  equation  (B16),  the  cosine  of  the  phase  difference  of 
the  equivalent  realization  is  equal  to  the  mean  of  the  cosines  of  the  phase 
differences  of  the  two  realizations  which  are  processed.  This  checks  with 
the  example  for  in  that  example 


A-  $ = 0° 

Bp  ap 


*Bq  - % “ 9°° 


then  from  equation  (B16)  we  have 


cos(<(>  - ) - H 

ar  Br 


so 


Wo  - $ ) - 60° 

Br  Yar 


(B23) 


(B24) 


